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Abstract 

We present a detailed analysis of the use of heavy quark fragmentation into heavy hadrons 
for testing the heavy quark effective theory through comparison of the measured fragmentation 
parameters of the c and b quarks. Our analysis is entirely model independent. We interpret 
the known perturbative evolution in a way useful for exploiting heavy quark symmetry at low 
energy. We first show consistency with perturbative QCD scaling for measurements done solely 
with c quarks. We then apply the perturbative analysis and the heavy quark expansion to 
relate measurements from ARGUS and LEP. We place bounds on a nonperturbative quark mass 
suppressed parameter, and compare the values for the b and c quarks. We find consistency with 
the heavy quark expansion but fairly sizable QCD uncertainties. We also suggest that one might 
reduce the systematic uncertainty in the result by not extrapolating to low z. 
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1 Introduction 



Current e + e~ colliders have excelled in the precision study of weak interactions. One can also hope 
to perform detailed tests of the predictions of QCD and heavy quark symmetry. With the copious 
production of heavy quarks, the study of heavy quark fragmentation functions at e + e~ colliders might 
prove one such forum for detailed studies of QCD. 

The factorization theorem states that the measured fragmentation function can be written as a 
convolution of two terms, one which depends on short- distance physics and one which is sensitive to 
large-distance dynamics. The perturbative contribution has been understood for some time, while 
heavy quark symmetry provides constraints on the nonperturbative contribution. 

The nonperturbative analysis that we follow was derived in ref. [|I|. The basic idea is to evaluate 
the matrix elements of the heavy quark operators at sufficiently low renormalization mass scales, 
where one can exploit heavy quark symmetry to expand the moments of the fragmentation function 
in A/m, where A is some QCD related scale, m is the heavy quark mass, and the coefficients in 
the expansion are mass independent. With this parameterization of the nonperturbative features of 
the fragmentation function at low mass scale, one can evolve the moments via perturbative QCD 
to the high scale Q 2 at which they are measured. One can then use the experimental results to 
extract the unknown nonperturbative parameters. Since the same parameters enter the c and b quark 
fragmentation functions and moments, one can test how well the heavy quark expansion converges. 

In this paper we apply this approach to the normalized second moments of the c and b quark 
fragmentation functions, measured at different center of mass energies. We study the second moment 
because it is best measured experimentally and because it involves a single nonperturbative parameter 
at order A/m. 

Our analysis consists of two parts. We first relate the second moments of the c quark fragmentation 
function measured at different center of mass energies, namely by ARGUS or CLEO and at PEP and 
LEP. This part of the analysis involves solely perturbative QCD evolution. In the second part of the 
analysis we combine both perturbative QCD and the leading order heavy quark expansion to relate 
the moments of the c and b quark fragmentation functions. 

We will see that there is no inconsistency with preliminary measurements. We will also see that 
by comparing the fragmentation function of c quarks at different center of mass energy one might 
hope to constrain the QCD scale. When using heavy quark symmetry to relate the b and c quark 
fragmentation functions, the calculation is less predictive, due to the large uncertainty in higher order 
QCD effects at the heavy quark mass scale. How useful the relations ultimately prove is sensitive to 
the size of Aq C d- 

In Section 2 we describe in detail the perturbative QCD scaling accurate to subleading logarithmic 
order, and we summarize the nonperturbative results of ref. relevant for the purposes of this paper. 
In Section 3 we present our results. Section 4 is a digression from our main analysis in which we 
suggest that the experimental results be presented differently. Rather than extrapolating the data to 
low z, we suggest a "cutoff" moment, where one only integrates from some cutoff value z . We suggest 
that factorization, though no longer exact, might work adequately. We conclude in the final section. 
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2 Theoretical Background 



In this section, we give the theoretical framework with which we analyze heavy quark fragmentation. 
We employ perturbative QCD in conjunction with the heavy quark expansion to relate the heavy quark 
fragmentation functions measured at different center of mass energies. Because it is best measured and 
involves only a single nonperturbative parameter (if we only keep the leading quark mass dependence), 
we will ultimately relate not the fragmentation function itself, but the normalized second moment. 

The perturbative analysis uses refs. f2|, || and the nonperturbative analysis is from ref. 0. From 
the latter reference, we use only the results that the moments of the fragmentation function have a 
heavy quark mass expansion, and that the correct choice of fragmentation variable corresponds to the 
ratio of energies. The discussion which follows primarily summarizes the perturbative scaling which 
was not done explicitly in ref. We frame the analysis entirely in terms of the standard procedure 
which is used when applying the heavy quark theory, namely scaling, matching to the heavy quark 
theory, and evaluating matrix elements in the heavy quark theory. With the perturbative scaling, the 
measured fragmentation function will be used to determine the unknown nonperturbative parameters 
of the fragmentation function. The tests of the theory come when we relate the fragmentation function 
measured at different center of mass energies, and when we relate the fragmentation parameters for 
different flavors of heavy quark. 

In this paper, we work to subleading logarithmic order. We do not incorporate explicitly Sudakov 
logarithms in the nonperturbative analysis, since they are already included in the parameterization of 
the heavy quark fragmentation function at low renormalization mass scales. Furthermore, the heavy 
quark expansion constrains us to the study of low moments, because of the large combinatoric factor 
which multiplies terms in the expansion of higher order moments. Therefore, as argued in ref. |J, 
Sudakov logs are not important in our approach, since low moments are not dominated by z in the 
regime where Sudakov logs are large. 

We use the following notation: a is the cross section for the inclusive e + e~ annihilation into a 
hadron, H, <7j is the short-distance cross section for producing a parton i, and fa is the fragmentation 
function for producing the hadron H from the parton i. According to the factorization theorem, the 
cross section for the inclusive annihilation, e + e~ — > H + X, can be written as 

da(z,Q,m) ^ f 1 dy d&i(y,Q, (i) ~ 

- A = 2W j fi{z/y,fi,m) (1) 

dz i Jz V dy 

where z is the ratio of the energy of the hadron to the beam energy in the center of mass frame 
(0 < z < 1). The factorization formula ([]]) states that the physical cross section is the convolution 
of two terms. The first, d&i(y , Q , //) / dy , is the short-distance cross section for producing a parton i 
(where % can be a heavy or light quark or antiquark, or a gluon). It is insensitive to the low-energy 
features of the process and therefore does not depend upon the mass of the heavy quark. 

The inclusive short-distance cross section for producing a parton % in e + e~ annihilation via 7- and 
Z-exchange, is given by 



dai(z,6,Q) 
dzd cos 9 




where z = 2q ■ p/q 2 , p^ is the four-momentum of the parton, is the four-momentum of the virtual 
boson, with y/q 1 = Q, and 9 is the cm. angle between the parton and the beam. We have taken fi = Q, 
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so that da q (z,9,Q) is /i-independent. The notation Ci,i(z) = Ci } r{z) and Ci^(z) = Ci y r{z) — 2Cj ) t(z) 
is also used in the literature. 

In eq. (0), the factor X q for quarks is just the total cross section for producing a quark 

K = ^ = -qT 1 ^ + 2e e e g v e v q Re X (Q 2 ) + (v 2 e + a 2 e )(v 2 + a 2 q )\ X (Q 2 )\ 2 } (3) 

where 

X(g2) = Q>-M*+iQ*T z /M z (4) 

and e e , u e , a e , e q , v q , a q are the charges and Z couplings of the electron and quarks, respectively. For 
gluons, Xq is given by 

Ac = + * ? ) (5) 

Integrating (0) over the angular variables we obtain 

L^i = Ci (z) = C i>T (z) + C i}L (z) (6) 
Aj az 

The coefficient functions Cj are calculable in perturbative QCD as a power series in a (//.): 

Ci(*, Q, /i) = Cf\z) + ^-Cf\z, Q,fi) + ... (7) 

The leading order term in the quark coefficient function C^l{z) is suppressed by a factor m 2 /Q 2 and 
therefore negligible at high energies, but the 0(a) contribution is not suppressed, so at subleading 
order we need to include both contributions to C q . The coefficient function for gluons, Cq, is of order 
a. 

The second factor appearing in equation (|1|) is the fragmentation function fi(z, fi,m), which de- 
scribes how the parton i combines with surrounding partons to produce the observed hadron H, in 
a much longer lapse of time. It is insensitive to the high-energy part of the process and therefore 
depends on m but not on Q. 

The scale \i is an arbitrary scale which separates the low- from the high-energy dynamics. If we 
know fi(z,fio,m) at some scale /io, we can obtain fi(z,fi,m) by using the Altarelli-Parisi evolution 
equations, which organize the large logarithms log(/x//xo) to the desired order in perturbation theory. 
In practice, we will work to subleading order, using the calculated anomalous dimensions from ref. 
0. As we will see, Ci(z,Q,/jb) contains terms of the form logQ 2 //i 2 , so we will choose /x 2 ~ Q 2 to 
avoid large logarithms. On the other hand, we take /io of order of the heavy quark mass, m, where 
the results of the heavy quark effective theory apply. 

It is convenient to define the singlet (/ ) and non-singlet (j^ 5 ) linear combinations of the quark 
fragmentation functions: 



(8) 

— f S 
2N P J 



H+) 

Jq 


= f. + fq 




N F 




= E/i +) 




9=1 


fNS 
Jq 


_ f(+) 
J q 
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where Np is the number of flavours and / denotes the fragmentation function of the antiquark q. 

According to the result of ref. 0, we can expand /q(z, fi Q ,m), the heavy quark fragmentation 
function, for // « m, Q as 



m \ /I m 



/ Q (z, m, m) = -a \^-Mil ] + $ [ | (9) 

where M# is the mass of the heavy hadron and e ~ 1 — -ff- . The funny dependence on z is because the 
expansion is simplest in terms of the variable 1/z. It is important to use z = E/Ef, eam (or its inverse) 
as the fragmentation parameter. Other fragmentation variables will permit nonvanishing higher twist 
contributions at leading order in the heavy quark mass expansion and oiqcd [01- With this choice of 
fragmentation variable, higher twist effects should be negligible (this is of course only significant for 
relatively low Q 2 (relative to m 2 ) experiments, like ARGUS and CLEO). 

We will not use the full form of this result, as the measured fragmentation functions are too poor 
to do a detailed fit, which involves many nonperturbative parameters. Instead, as suggested in ref. 
|TJ], we will consider only the moments, defined as 

V m) = f dzz N - l f Q {z, /i, m) (10) 

J 

The low order moments have the advantage of being better measured, and of being describable in terms 
of a small number of nonperturbative parameters (when the heavy quark expansion is exploited). 

According to the assumptions of the heavy quark theory, the fragmentation functions(and their 
moments) for the light quarks and gluons at the low renormalization scale can be taken to vanish. 

By taking the moments of the function of eq. @, we derive a heavy quark expansion for each of 
the moments. In general, at order A/m, there are two independent parameters in terms of which any 
moment can be expanded with known N dependence [[!]. However, since we will restrict our numerical 
analysis to the measured second moment, it is convenient to define a single nonperturbative parameter 
and expand the second moment as 

ff Q (m,m) = l-- (11) 
m 

From previous work 0, one can also conclude that there is a heavy quark mass expansion for the 
moments. The power of the approach in [Q] is that one can address the higher twist corrections to 
determine the best choice of fragmentation variable, that one can relate higher order moments, and 
the one can in principle include higher order mass suppressed terms. The parameter a is the single 
nonperturbative parameter which is required to describe the second moment of the fragmentation 
function if we work at order A/m. Notice that because the fragmentation function is not defined in 
terms of an on shell matrix element, the parameter a is not determined theoretically, but can be fit. 
(This is in contrast with the distribution function, for which one can derive a(x) = 5(x — m/M), with 
corrections of order (A/m) 2 .) In ref. [pj, a was defined in units of A = M# — m. In this paper, we 
have absorbed A in the definition of a, which is now a dimensional parameter. It should be kept in 
mind however, that a is expected to be of order Aq C d- The parameter a depends on the hadron type. 
Higher order terms can be incorporated-however if the heavy quark expansion is valid, these terms 

1 This expansion is true in the HQET, so Jq denned here is related by a nontrivial perturbative matching condition 
to /q used in eq. (|l|). 
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should be suppressed, and we should get an answer accurate at about the 10 % level by just keep the 
first mass suppressed contribution. 

This expansion is valid only for /xq ~ m. It is therefore necessary to match onto the full theory, 
and scale the result to \i ~ Q. We now describe the perturbative scaling and matching procedure. 

The fragmentation functions obey the Altarelli-Parisi evolution equations 

dffi^ 2 m) = Ef^/^))rW,™)- tS = s,g (12) 

a log// 2 ^Jz y 

df" s (z,fi,m) [i m dy 



., P N8 (z/y,a(jj))f^(v,fi,m)-^ (13) 

alog/i z Jz y 

Only singlet quarks and gluons mix under evolution []; the non-singlet equation ( [□!) decouples and 
can be solved independently. 

The evolution functions P TT '(z, £*(//)) and Pns( z , «(/-0) can be expanded in a power series in a(fi): 

P(z,afr)) = ^-PW(z) + (^) 2 P {1 \z) + 0(a 3 ) (14) 

The leading log contributions {P^'i Pns = ^ss) are the well-known Altarelli-Parisi splitting functions 
1- 

The definitions of Cf and fi, and consequently of P(z, a (//)), are unique only at leading log order. 
At the next-to-leading level they depend on the factorization scheme used to separate the collinear 
singularities. There are complete calculations available in the literature in two different factorization 
schemes, which we will denote I and II for definiteness. Scheme I is defined in pi by Curci, Furmanski 
and Petronzio, where they give the quark coefficient functions and non-singlet evolution function in 
z-space, as well as their Mellin transforms to A-space. The evolution function matrix in z-space for the 
singlet sector, including quark and gluon mixing, has been computed in |7j and the gluon coefficient 
functions can be found in 0. Matching conditions for quarks and gluons within this scheme are 
calculated in ||. Alternatively, in || Floratos, Ross and Sachrajda define Scheme II, which is used in 
to compute coefficient functions for quarks and gluons as well as non-singlet and singlet evolution 
functions both in z-space and A-space []. 

The results for P^, in [0 and || do not agree; there is a difference for Pqq(z) in the term propor- 
tional to the colour factor Cq. The timelike evolution functions (14) are obtained from the spacelike 



ones (adequate to deep inelastic scattering processes) by analytic continuation. The discrepancy men- 
tioned above is in fact a consequence of the discrepancy in the corresponding spacelike evolution 
functions previously computed 0. 

In principle one should evolve the fragmentation functions fi(z,n,rn) by using eqs. (|i~2| ) and (|i~3|), 
i.e., including mixing between singlet quarks and gluons. However we have explicitly verified using 
Scheme I that incorporating mixing never modifies our result by more than about 1% for the b quark 
and 8% for the c quark. Therefore, for the purpose of an expansion to order A/m, we could safely 
neglect mixing. The reason mixing is small is that the gluon fragmentation function is never very 

2 Note that the anomalous dimension matrix for timelike processes is transposed aside from the dependence on Nt. 
3 However, expressions for the timelike moments of the coefficient functions and anomalous dimensions given in the 
appendix B of ref . || correspond to and , despite the N + f index. 
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big and the gluon coefficient function is O(a). Our result agrees with the fact that the probability of 
producing a heavy quark indirectly, via a secondary gluon, is very small at Q 2 ~ M|. The authors 
of ref. [I(J have shown that the so-called Branco et al. terms, which contribute to the heavy quark 
inclusive cross section, can not be larger than 0.5% for 6-quarks and 4% for c-quarks and one would 
expect that the remaining terms are of the same order of magnitude. This is very important to our 
analysis, since it means that we do not need to scale the first moment (or equivalently the normalization 
factor) between different momentum scales. If we are only dealing with energies less than or equal to 
the Z mass, the normalization is approximately independent of energy scale; that is very few additional 
heavy quarks are produced through QCD processes. Therefore, we consider only non-singlet evolution 
in the remainder of this paper. 

Beyond leading log level, the non-singlet evolution function has an additional contribution due 
to qq mixing. However, we have found that it is negligible at the energies which are currently of 
interest. Therefore, we use the non-singlet anomalous dimensions given in ||, which correspond to the 
factorization scheme I and do not include qq mixing. 

We are interested in the evolution of the moments of the heavy quark fragmentation function, 
Tjv(/i, m). The next-to- leading order evolution equation for Tjy has a very simple form: 



df N (n,m) 
d log ji 2 



P N {a{ix))t N {ix,m) 



(15) 



where 



P N (a(ii)) 



a (v) p(0) , 



'«(//) 
2tt 



P 



(i) 



N 



(16) 



is the Mellin transform of the Altarelli-Parisi splitting function in next-to leading order. The formula 
for Pff is given in the appendix of ref. ||. The expression for a(fi) accurate to next-to-leading order 
is 

1 / 6 1 log[log(/i 2 /A 2 )]\ 



with 



Then eq. ([15]) takes the form 

dT N (fji,m) 



d log fi 2 



b log(/i 2 /A 2 ) 

33 - 2N F 
" 12k 



P (o) , 



2tt 



b 2 log(/i 2 /A 2 ) 

153 - 19iV F 
24^ 



;i7) 



p(l) 27T&! (o) 

N ~bV N , 



t N (n,m) 



(19) 



which can be solved analytically, yielding 



t N (fi,m) = T N (no,m) 



Oi(fio) 



p(0) 

fjv 



exp 



a(/io) - a(/x) 
4tt 2 6 



'p(i) H (o)' 



(20) 



The physical fragmentation function (eq. ([!])) is obtained by convolution of the fragmentation 
function with the hard cross-section, evaluated at the scale /i. Therefore, the normalized moments 



a N (Q,m) 



-f 

a Jo 



dzz N-i d(T ( z 'Qi m ) 



dz 



(21) 



6 



are given by 



,(») 



a N (Q,m) = C N {Q,n) 





-A— r 

2nb f 


exp < 





2 ^o f oQ/q) - a(/x) ^(i) 2tt6i 



47T 2 6 n 



3 (0) 
N 



where Cn(Q,h) is the Mellin transform of the coefficient functions (eq. (|7|)) 

C N (Q,fi) = 1 + ^C%\Q,n) + 0(a 2 ) 



2tt 



t N (no,m) (22) 



(23) 



The expression for C^\Q,fi) is 

^ 1} (Q,/i) 




c f 



+Sf(N) + S l (N) - 



2Si(JV) 



n(jv + i; 



5S 2 (iV) 



1 2iV + 1 1 

+--2 ^ n/ _ — + 



3 1 
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N N 2 (N + 1) 2 (N + l) 2 2iV + l 2 



where Cp = 4/3 and 



A? 



Si(N) = E 



(24) 



(25) 



3=1 



We will take \i 

(i) 



Q from now on, so that C^ 1 does not depend upon /x. We therefore use the notation 



C%'{n,n). 



In eq. (|2~2|), r^(fj, , m) is the moment of the heavy quark fragmentation function in the full theory 
at the scale //o ~ Jn. Note that it is critical to the use of a subleading log calculation that we know the 
/i-dependence of the matrix element. Even though the matrix element itself is nonperturbative, the 
scale dependence is not, and is contained in the matching coefficient. We note that in reality, there 
is still some scaling which goes on below the heavy quark mass scale because of the scale dependence 
of the spin-dependent mass suppressed operator []. However, since the maximum scaling between the 
b and c mass scales of the gluon magnetic moment operator is 20%, and it is only one of several 
contributions to the A/m term, this effect is never numerically more significant than the higher order 
mass suppressed corrections. We therefore neglect scaling in the heavy quark theory itself; the /i 
dependence is then just that from the matching between the full and heavy quark theories. Because 
this matching is done at a fairly high momentum scale, of order the heavy quark mass, the scale 
dependence is calculable order by order in perturbation theory. We then have 



r N (ii ,m) = 1 + 



«(^o) 
2tt 



d {1) 

a N 



(// ,m)J T^ Q (m) 



(26) 



We take f^(m) as independent of fio, but because we work at a finite order in perturbative theory 
we will see that there is a fairly substantial dependence on /j,q which cannot be neglected. This will be 
clear in section 0. 



4 C (D, 



N '{Q,n) = d$(Q,[i) in the notation of rcf. 
5 We thank Eric Braaten for pointing this out. 



7 



The matching condition which agrees with the factorization prescription for Cjp(Q,fi) and 



can be obtained from ref. p and reads 
dx\no,m) 



-2S 2 (N) - 



-2S*{N)+2S l {N) 1 + 



+ 2 



N(N+ 1) 



where 



(0) 



TV 



(iV + 1) 2 N(N + 1) 

2S X (N) 



r r i 1 

F ^2 iV(iV + l) 



Using (|26| ) we can write eq. (|22| ) as 



(T N (Q,m) = {z N -\Q,m)) = {z{Q,my-% ert {z{m) N - 1 ) 



N-l\ 



\N-1\ 



nonpert 



(27) 



(2f 



(29) 



where 



(z(Q,m) N ^pert 



(^(TTi) }nonpert 



C N (Q,/i) 



p (0) 
2ir6 () 



x exp -j 

-pHQ 

— TV ' 



a(Q) 
a(/i ) - a(Q) 



'1 , a (^o) ,(1)/ \ 
1 + -^-d K N J (tiQ,m) 



N 



2tt6i (o) 
bo 



m) 



(30) 
(31) 



It is important to clarify the issue of the factorization scheme dependence. The separate results for 
the anomalous dimensions and the coefficient functions depend on the factorization scheme, but in the 
convolution there is a cancellation of the scheme dependence order- by-order in a 0. Neglecting next- 
to- next-to-leading terms and using eq. ( p3f) , we can write the coefficient of a(Q) in the exponential of 
eq. (|30D as ~ 2irboCff . This combination can be shown to be independent of the renormalization 
scheme at the subleading level, and we have explicitly checked that by using both schemes I and 
II. However the coefficient of a(/io) still depends on the factorization scheme and therefore so does 
the matching condition (|27|). In order to obtain an expression in terms of renormalization scheme 
independent quantities, we use the equation 



(/, ,m) = i + ^4W) + ^4 1 > 

Z7T ZTT 



(32) 



where C^ Q (/io, m) is the short- distance cross section for producing a heavy quark of mass m. We then 
rewrite eq. (j30"D as 



(z(Q,m) N l ) peTt = ('" Q (/i { 
x exp 



Q(/ip) 
a(Q) 
a{n ) - a(Q) 



p(0) 

N 
2jrf)Q 



An 2 b n 



P 



(i) 



N 



2nb C 



(i) 

N 



27T&! 



P 



(0) 



N 



(33) 



8 



Eqs. (^) and are equivalent at the subleading level. Moreover, the evolution factor entering in 
eq. (|33D is now independent of the factorization scheme and thus the matching condition (|32|) , derived 
in the scheme of 0, is also adequate for the scheme used in ||. 

However, since we are using results calculated within a single scheme it is not necessary to use 
scheme independent expressions as in (^). In fact, this turns out to be the preferred procedure. This 
is because in order to obtain a scheme independent result, we added terms which are technically higher 
order, but are in fact quite large, due to the large size of C$. For example, ^f^C^ is as big as 
0.8 for A5 = 225 MeV. So we are probably introducing spurious /1 dependence by adding these unduly 
large terms. We therefore used eq. ( |30D in our calculations. 

To summarize, we have shown that the measured moment of the fragmentation function can be 
expressed as in eq. (|9|). The perturbative contribution (z(Q, m)) pert is given by (|30|) or (|33|), where the 
former has smaller /1 dependence and is actually what one obtains on integrating the renormalization 
group equation. Once one has calculated {z(Q, m)) pert , one can extract the nonperturbative parameter 
a, defined in eq. ([TTJ) . This parameter is independent of Q, the energy of the experiment, and so long 
as one always deals with the same kind of hadron, it should be also independent of m. 

Unfortunately however there are still fairly large QCD uncertainties in the perturbative scaling 
given that we are working only to subleading logarithmic order and that the QCD scale is not so well 
known. The problem is that we need to know (z(Q,m)) pert sufficiently well to extract (z(m)) nonper t 
at the level of a/m. For large m, the required accuracy is greater, though the perturbation theory 
is better. We find that the error in extracting a from the b and c fragmentation functions is in fact 
comparable. In order to determine the accuracy of our procedure, we will consider different values of 
jiQ. Higher order perturbative effects may also be estimated from the 0(a 2 ) terms neglected in eq. 
(P3|). We find that these uncertainties are never as large as the /zo dependence. So we estimate the 
uncertainty from higher order effects by considering different values for the renormalization scale /i . 
The particular choices and results are given in the next section. 



3 Results 

In this section, we apply the procedure described in the previous section, where it was shown how 
perturbative QCD and the heavy quark expansion can be used to relate moments of fragmentation 
functions of different flavors measured at different Q 2 values. These predictions can be compared 
with the data from ARGUS, CLEO, PEP, PETRA and LEP. Since current results from LEP are 
only preliminary, we present our analysis in such a way that it can be readily applied with improved 
measurements. To do this we center the values on our plots on the preliminary measurement and we 
extract the nonperturbative parameter over the full range of experimentally allowed numbers (that is, 
the measured value with 2 sigma errors). 

Recall the uncertainties in our predictions arise from two different types of QCD uncertainties. 
There is uncertainty in the perturbative part of the calculation due to both the poorly known value 
of oiqcd an d the fact that we work only to subleading logarithmic order. For quantities involving 
perturbative scaling only for large Q 2 these uncertainties are expected to be small. But in order to 
apply heavy quark symmetries, one always needs to scale to the quark mass scales. For such low energy, 
higher order QCD corrections and the uncertainty in the exact value of oiqcd can be important. 

It therefore makes sense to divide our analysis into two parts. We first consider simply measuring 
(z c ) at two different energy scales (sufficiently larger than the c quark mass that higher twist contribu- 
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tions should be small). These values should be related by only perturbative QCD, at scales for which 
the subleading calculation should prove reliable. The measurements seem to be consistent with per- 
turbative QCD scaling over the allowed range of Aqcd- Turned around, it means that with sufficiently 
accurate data, one can constrain Aqcd (although probably not as well as with other methods). 

Following this analysis, we proceed to incorporate the full formalism, employing both perturbative 
QCD and a nonperturbative expansion to relate b and c quark fragmentation functions. Here we will 
find sizable QCD uncertainty. Nonetheless, we will see that the data is substantially self consistent. 

Notice that throughout the analysis of section |2| we have assumed fragmentation into a specific 
final state, although experiments at very high center of mass energy such as LEP do not distinguish 
hadron species. This should not be a problem when comparing inclusive measurements such as those 
done at LEP, since the normalized second moment for the inclusive measurement has the same heavy 
quark expansion for the b and c quark. 

This is potentially a problem if different experiments select differently on the final hadronic state. 
The problem occurs because the nonperturbative parameter a depends on the hadron type. The 
nonperturbative contribution to the mean (z) experimentally measured can be a different linear com- 
bination of the different a's at the different experiments. However, measurements from ARGUS show f] 
that (zd) ~ (zd*) ~ ( z a c ) within errors. Preliminary measurements from CLEO seem to also support 
this claim. If this is the case, we can to a good approximation neglect the differences between the 
different hadron types when comparing the measurements at different center of mass energy. 

Moreover, for inclusive measurements at sufficiently high energy, the dependence on hadron type 
can be ignored, since the perturbative QCD scaling is the same for all hadrons so that it is always 
the same linear combination of a's which is measured at any energy. This is obvious, since /i can be 
chosen so that only the coefficient function depends on Q. 

Finally, for measurements which select on a specific final state (eg D*), there is no problem. 

Before we begin, we outline the differences in our approach from previous studies [O. First, we 
make no attempt to fit the entire curve. We look only at the mean (z), which is much better measured 
than the fragmentation function itself. Furthermore, the mean (z) depends only on QCD parameters 
and one single nonperturbative parameter a, where a is defined by 



Our result is therefore independent of any assumed functional form of the fragmentation function. 
Furthermore, (z) nonper t is much more stable against variations in \x and the QCD scale than the 
parameters used to fit a particular functional form. For example, in ref. JTTJ, where the functional 
form z nonpert — z a (l — z) 13 was used, the parameters a and (3 varied enormously in comparison with 



Second, we use eq. (|30|) rather than (|33|) to calculate the perturbative factor (z(Q,m)) pert . As 
explained earlier, this result is better behaved with respect to /i-variation, and it is in fact what is 
obtained by integrating the renormalization group equations. 

Third, we do not incorporate explicitly Sudakov logs in the nonperturbative fragmentation function. 

Fourth, we use the experimental results on the mean (z) = (E) / Eb eam . It is important to use the 
variable (z) because its non-perturbative contribution scale linearly in the mass of the heavy quark, 
according to eq. fl3~4]). 

6 We thank G. Bonvicini for analyzing the data from the Ph.D. Thesis of P. C-Ho Kim and H.C.J. Seywerd (ARGUS 
Collaboration) . 




(34) 
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Finally, we use updated experimental data from LEP which leads to a harder b quark fragmentation 
function. As a result we find better agreement with the heavy quark effective theory predictions than 



It should be noted that there is some dependence in the measured (z) value on the assumed 
functional form (contributing to the systematic error) because the data is not measured to arbitrarily 
low z, so functional forms of the fragmentation function are assumed in order to extrapolate to low 
z. In section |] we will show that one can use moments with a cutoff greater than zero to reduce this 
dependence on the functional form. 

In both parts of the following analysis, we will allow for different values of Aqcd within the allowed 
range. We estimate the importance of the neglected higher order perturbative QCD contributions by 
testing the stability of our results with respect to a change in renormalization scale. It is an important 
feature of the analysis that the /io dependence is entirely perturbative for a sufficiently heavy quark. 
Nonetheless, a(/io) at the scale of the heavy quark mass is so large that it is not clear how well the 
perturbation theory converges. We estimate this uncertainty by allowing for a sizable variation in /i , 
between 2m and m/2. This might be too large a range, but without a higher order calculation, it is 
impossible to determine the accuracy of the subleading calculation. 

When we invoke the heavy quark expansion, we work in the approximation of only retaining 
the leading mass dependent correction, suppressed by a single power of the heavy quark mass. Of 
course one can readily incorporate higher order terms. But we would like to test whether the heavy 
quark expansion works sufficiently well to get agreement (at least at the 10 % level) even without 
incorporating higher order contributions. This is important, as the heavy quark effective theory is still 
essentially untested at the A/m level. 

3.1 Perturbative Analysis: Relating (z c ) at Different Energy Scales 

In this section we focus on what can be learned using solely perturbative QCD evolution, without 
the heavy quark effective theory. Even if it turns out that QCD perturbation theory at subleading 
order is not sufficiently accurate at the low scale /Zq ~ m 2 , it certainly should give a good description 
of the evolution between scales above the mass of the T. Thus if we compare the moments of the 
c quark fragmentation function measured by CLEO or ARGUS, and at PEP, PETRA and LEP, 
the measurements should be related through the perturbative evolution between the scales of the 
experiments. Again, this is with the caveat that if a different linear combination of hadron states is 
selected, the relations only need hold if the nonperturbative fragmentation function is the same for 
each hadron type, as seems to be the case. Furthermore, since D* J s are the dominant contribution 
at these experiments, one would expect the perturbative relation to hold fairly well. In fact, with an 
exclusive measurement on D*'s at both experiments, the perturbative scaling result would certainly 
apply. 

According to the evolution equations, the N = 2 measured moments of the fragmentation function 
for the c quark are given by 



in ref . ITTt . 




(35) 



Therefore the ratio between the mean (z c (Q)) measured at different energy scales is 



(z c {Q 2 )) _ C 2 (Q 2 ,fi) 
(z c {Qi)) C 2 (Q uf j) 



(36) 
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A / 1\ IT T T \ 

A 5 (MeV) 


H = Qi/2 


H = Qi 


H = 2Qt 


75 


0.813 


0.825 


0.832 


125 


0.794 


0.808 


0.816 


175 


0.779 


0.795 


0.804 


225 


0.766 


0.784 


0.794 


275 


0.755 


0.774 


0.785 



Table 1: = 91GeV,„) 



C 2 (Q 1 = 10.55GeV,/i) 

We relate the N = 2 moments measured at ARGUS (Qi = 10.6 GeV) or CLEO (Q 1 = 10.55 GeV) 
to those measured at PEP (Q 2 = 29 GeV) and LEP (Q 2 = 91 GeV). Since the experimental data 
from ARGUS and LEP have smaller errors, we focus our discussion mainly on the results from those 
experiments. The difference between the perturbative calculation for ARGUS and CLEO energies is 
negligible, so we only present the theoretical results for Qi = 10.55 GeV. 

Since for the theoretical computation of the right hand side in eq. (^) we only need to evolve 
the coefficient functions between Q\ and Q 2 , which are well above Aqcd, we expect the perturbative 
evolution with next-to-leading accuracy to work well. The evolution equations for the coefficient 
functions are related to that for the fragmentation function by 

fj J dC N (Q,fj J ) _ ndf N (fi,m) 
dfi dfi 

We choose the scale /i ~ Qi, so C 2 (Qi,A*) is given by eqs. ([2"3]) and (^|), with Q = Q 1 . For 
^2(^2,/^) we also use the evolution equation fl22| ) and eq. ([37]) to obtain 

,(0) 



a(/x) 
a(Q 2 ) 



To account for theoretical uncertainties of the prediction, we computed the ratio C 2 (Q 2 , fi)/C 2 (Qi, ji) 
using the values of A in the five flavor theory, A 5 = 75, 125, 175, 225, 275 MeV and, as an estimate of 
higher order effects, we varied the scale \i between Qi/2 and 2Qi. 

The QCD perturbative results for the ratio C 2 (Q 2 , /i)/C 2 (Qi, fi), with Q 1 = 10.55 GeV and Q 2 = 91 
GeV are shown in Table 1. Notice that the results do not depend strongly on the scale /1, as one would 
expect since subsubleading effects at the energy scales we are considering should be rather small. 

The mean (zjj*+) for the fragmentation function of a c quark into D* + mesons has been measured 



both at ARGUS 12 and LEP ff3|. The results are 



(z D *+{Qi = 10.6GeV)) = 0.647 ±0.006 
(z D *+(Q 2 = 91GeV)) = 0.495 ±0.013 (39) 

In order to relate the inclusive measurements from PEP and LEP with the ARGUS measurements 
of (zd), (zd*) and (za c ) we have to average over all the final states. [] If we assume heavy quark 

7 In fact, there are ARGUS measurements only for D* + mesons, but we assume (z^.o) = (z D *+). 
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Figure 1: (z c (Q2 = 91GeV)) as a function of (z c (Qi = 10.55GeV)), for A 5 = 125 MeV (dashed line), 
175 MeV (solid), and 225 MeV (dashed-dotted). For each value of A5 the lower line corresponds to /i — Q\/2 
and the upper line to fl — 2Q\. 



symmetry and that a c quark fragments into baryons about 12% of the time, the mean value of z c is 



[z c {Q x = 10.6GeV)) 



0.12(z Ac ) + O.i 
0.640 ± 0.009 



-(z D ) + — (zd* ) 



(40) 



It would be very useful to test the assumed ratio of 3:1 for D* relative to D production is in fact 
correct. Preliminary results from CLEO seem to give a bigger value, so we allow for a wide range 

is 



0.62 < (z c (Qi)) < 0.68. The mean (z c ) measured inclusively at PEP (14 



(z c (Q 2 = 29GeV)) = 0.526 ± 0.03 



(41) 



while at LEP 

{z c (Q 2 = 91GeV)) = 0.487 ± 0.011 (42) 

The numbers shown in Table 1 are all that is needed to compare the measurements from ARGUS or 
CLEO and LEP. In Figure 1 we present this result graphically for a slightly narrower range of A5, 
namely A 5 = (175 ± 50) MeV, according to ref. fli~5| . For each value of (z c (Qi = 10.55GeV)) given 
on the horizontal axis, we plot the value of (z c (Q2 = 91GeV)) which should be measured according to 
eq. ([J6|). There are two lines for each value of A 5 , corresponding to \i = Qi/2 and fi = 1Q\. So for a 
given value of cxqcd the prediction for (z c (Q 2 = 91GeV)) lies between these lines. 

We see that from the ARGUS measurement (z c {Qi = 10.6GeV)) =0.640, (z c (Q 2 = 91GeV)) is 
predicted to be in the range 0.490 < (z c (Q 2 = 91GeV)) < 0.522. If we allow (z c {Q x = 10.6GeV)) to 
vary within 2a around the measured value we obtain 0.477 < (z c (Q 2 = 91GeV)) < 0.537. Therefore, 
the results are consistent at the 2a level but they seem to favour large values of Aq C d- 

If we restrict the analysis to D* + mesons, from (zd*+(Qi = 10.6GeV)) = 0.647 we predict 0.496 < 
(z D ,+ {Q 2 = 91GeV)) < 0.528 (0.487 < (z D *+(Q 2 = 91GeV)) < 0.538 within 2a). It is very important 



13 



to relate the exclusive measurements on D* + mesons at both experiments, because in this case the 
nonperturbative fragmentation function is exactly the same (as it does not depend on Q 2 ) and we are 
sure to be testing perturbative QCD, without invoking heavy quark symmetry at all. With sufficient 
statistics, this could prove to be a useful measurement of Aqcd- 
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<z (10.55 GeV)> 
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Figure 2: {z c (Q 2 — 29GeV)) as a function of (^ c (Ql — 10.55GeV)), for the same values of A5 and fi as 
in Figure 1. 

We have also related (z c (Qi = 10.6GeV)) and (z c (Q 2 = 29GeV)). The result is shown in Figure 2, 
for the same values of A 5 and \i as in Figure 1. In this case we obtain 0.549 < (z c (Q 2 = 29GeV)) < 0.578 
from the ARGUS measurement (z c (Q 1 = 10.6GeV)) =0.640 (0.534 < (z c (Q 2 = 29GeV)) < 0.594 
within 2a). We also see that although the QCD perturbative calculation tends to give a too large 
prediction for (z c (Q 2 = 29GeV)), since the experimental error of this measurement is bigger than at 
LEP, the theoretical prediction lies always within 2a of the measured value. Therefore, it is not possible 
to constrain further Aq C d by relating (z c (Qi = 10.55GeV)) and (z c (Q 2 = 29GeV)). Obviously, the 
same applies when we compare (z c (Qi = 29GeV)) and (z c (Q 2 = 91GeV)). 

One can do a similar analysis for the b quark, relating {zb(Qi = 29GeV)) measured at PEP and 
(zb(Q 2 = 91GeV)) from LEP. However, the QCD perturbative evolution is in fact inconsistent with 
the experimental measurements 



rj, |TJ, namely 

(z 6 (29GeV)> = 
(z 6 (91GeV)> = 



0.715 ±0.03 
0.714 ±0.012 



(43) 
(44) 



Presumably, the b measurement at 29 GeV is not reliable. 

We conclude that perturbative QCD predictions are consistent with current data on c quark frag- 
mentation function within experimental errors. However the data seems to favor large Aqcd- Figure 1 
shows that a measurement of (z c ) from CLEO would be very useful in constraining further Aqcd- 

The measured second moment of the b quark fragmentation function at Q = 29GeV seem to be 
inconsistent with perturbative QCD and the LEP measurement at Q = 91GeV. For the c quark, the 
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A 5 (MeV) 


Ho = m c /2 


Ho = m c 


/i = 2m c 


75 


0.842 


0.865 


0.888 


125 


0.794 


0.832 


0.865 


175 


0.743 


0.800 


0.845 


225 


0.685 


0.769 


0.826 


275 


0.618 


0.736 


0.809 



Table 2: (z(Q,m)) pert for Q = 10.55 GeV, m = 1.5 GeV. 



second moment measured at Q — 29GeV does not provide further bounds, since the experimental 
error is larger than in later experiments. Therefore, in the following section we only consider the 
more accurate measurements of the heavy quark fragmentation functions at center of mass energies 
Q = 10.55 GeV and Q = 91 GeV. 



3.2 Relating and (z c ) 

In this section, we incorporate QCD directions and the leading order heavy quark expansion in order 
to relate the measured second moment of the b and c quark fragmentation functions, the first measured 
at LEP and the latter measured at ARGUS, CLEO and LEP. 

From Section 2, we know that for each value measured at a particular value of Q 2 , 

(z(Q,m)) = (z(Q,m)) pert (z(m)) nonpert . (45) 

We find the first term via the Altarelli-Parisi evolution we described in Section and use the measured 
value (actually a range of possible values) to determine the nonperturbative factor. 

In Tables 2-4, we present the factor (z(Q, m)) pert for Q = 10.55 GeV, m = 1.5 GeV; Q = 91 GeV, 
m = 1.5 GeV; and Q = 91 GeV, m = 4.5 GeV, where we have taken /i = 2m, m, m/2 and the values 
of A5 shown in the first column. With these perturbative factors, we can extract (z(m)) nonper t from 
measurements of the fragmentation function at CLEO and ARGUS and LEP (providing the subleading 
calculation is adequate). 



The measured values at ARGUS and LEP are the following fl2|l , [|I 



(z c (10.6GeV)) = 0.640 ±0.009 
(z c (9lGeV)) = 0.487 ±0.011 (46) 
(* 6 (91GeV)) = 0.714 ±0.012 

The numbers from LEP are still preliminary. Because the exact numbers are not yet known, and in 
order that our analysis can be applied when more exact numbers are measured, in each case we extract 
the parameter a over a range of values, namely 

0.62 < (^ c (10.55GeV)) < 0.68 

0.46 < (z c (9lGeV)} < 0.51 (47) 
0.69 < (z fe (91GeV)) < 0.74 
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Figure 3: Non-perturbative parameter d (in GeV) as a function of (z(QjUl)) for Q — 10.55GeV and 
m = 1.5GeV, with A 5 = 125 MeV (dashed line), 175 MeV (solid), and 225 MeV (dashed-dotted). For each 
value of A5 we have taken Ho = 2t71, ?7T, TTl/2, the higher line corresponding to the larger jJL§. 



A 5 (MeV) 


A*o = rn c /2 


Ho = m c 


Ho = 2m c 


75 


0.695 


0.713 


0.733 


125 


0.641 


0.672 


0.699 


175 


0.590 


0.636 


0.671 


225 


0.537 


0.603 


0.647 


275 


0.478 


0.570 


0.625 



Table 3: (z(Q,m)) pert for Q = 91 GeV, m = 1.5 GeV. 

Using the perturbative factors from Tables 2-4, in Figures 3-5, we present (on the vertical axis) 
the resulting values of a (in units of GeV) which would be extracted if the value on the horizontal 
axis is measured. We have presented the data graphically, rather than in tabular form to illustrate the 
sensitivity to the exact value which is measured. The values of Q 2 and m are as in the tables. 

In each figure we give the extracted parameters for A 5 = 125 MeV (dashed line), 175 MeV (solid 
line), and 225 MeV (dotted line) (the current range in the particle data book Hi~5|). For each value 
of A5 we present three lines, corresponding to Ho — 2m,m,m/2. The largest values of a in each case 
correspond to the largest value of Ho- This is readily understood, since there is more QCD scaling for 
smaller /i , so a smaller a is required to get a particular measured value. 

By way of illustration, we give the predictions for a corresponding to A 5 = 175 MeV. First consider 
the mean value of the measurement. In Figure 3, one can see that from the ARGUS measurement 
(z c (10.6GeV)) = 0.640 we predict 208MeV < a < 365MeV . Similarly, from the LEP measurements 
(z c (91GeV)) = 0.487 and (z 6 (91GeV)) = 0.714 in Figures 4-5 we obtain 264MeV < a < 413MeV 
and 425MeV < a < 552MeV , respectively. Therefore we find no overlap between the a parameter 
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Figure 4: a (in GeV) as a function of (z(Q, m)) for Q = 91GeV and m = 1.5GeV. The values of A 5 and 
Hq are as in Figure 3. 



A 5 (MeV) 


Ho = m b /2 


Ho = m b 


Ho = 2rn b 


75 


0.829 


0.836 


0.846 


125 


0.807 


0.816 


0.828 


175 


0.789 


0.800 


0.814 


225 


0.772 


0.786 


0.802 


275 


0.757 


0.773 


0.791 



Table 4: (z(Q,m)) pert for Q = 91 GeV, m = 4.5 GeV. 

determined from the measured mean values (z c ) and (z b ). For A 5 = 125 MeV there is no overlap 
either, and for A5 = 225 MeV there is overlap only between the LEP measurements. However, if 
we allow for a la variation around the (z) measured in each experiment there is substantial overlap: 
from the ARGUS measurement of (z c ) we predict 175MeV < a < 396MeV , and from (z c ) at LEP 
208MeV < a < 461MeV , while from (z b ) we obtain 287MeV < a < 686MeV . 

There are certain qualitative features of agreement which are good to note. First, the parameter a 
is never larger than 800 MeV, and is very likely smaller (especially if the value is indeed in the overlap 
region of the c and b quark results). This is reassuring, as one could not presume to do a heavy quark 
expansion for the c quark for a parameter a much larger than this. Furthermore, we see there is indeed 
substantial overlap for each possible value of Aq C d between the allowed range of a for each experiment. 
Remember, according to the heavy quark effective theory, one would predict the parameter a is the 
same in each case up to higher order corrections, of order (A/m) 2 , which one expects to be of order 
10%. It is also clear that the allowed values of a as determined from the two (z c ) measurements are in 
very good agreement. Both LEP measurements have greater overlap than the ARGUS measurement 
has with (z b ), therefore a measurement from CLEO would be very useful. 
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Figure 5: a (in GeV) as a function of {z(Q, TTl)) for Q = 91GeV and 777. = 4.5GeV. The values of A5 and 
/irj are as in Figure 3. 

However, the value of the nonperturbative parameter a which is extracted is clearly strongly de- 
pendent both on the precise value of Aqqd and the renormalization scale /xq, indicating the importance 
of subsubleading logarithmic corrections. This is because the difference between the measured value 
and the perturbative contribution is much more strongly /i -dependent (defined as the fraction of the 
value of the quantity) than the perturbative contribution itself. So for example, if (zi,) measured at 
LEP is 0.71, the parameter a for A 5 = 125 MeV, can vary between about 537 MeV and 637 MeV, 
even though the perturbative factor varies only between 0.81 and 0.83. In order to determine an effect 
of size Aqc,d/^, the perturbative contribution needs to be known at accuracy more precise than this. 
Clearly, the subleading calculation is not adequate for a precise extraction of a. This problem is not 
endemic to the particular test of the heavy quark theory we have done. It seems likely that in order 
to determine A/m corrections, one either must take ratios in which the QCD corrections cancel, or do 
an even higher order calculation. 

However, the parameter {z) nonpert is more stable against QCD uncertainties. This is actually the 
quantity upon which the measured value depends. When we relate two different measured values of 
(z), the variation with /x and Aqcd might not be so large. That is, we have not at all correlated the 
choice of fio/m for the b and c quarks. 

In Figures 6-7, for the assumed values of (z c (Q)) which are given on the horizontal axis, we 
determine the nonperturbative parameter a. We then use this to predict the value of (z&(91GeV)) 



which should be measured for the b quark (this was the procedure followed in ref. JTTJ] but for the 
fragmentation function itself). There are two graphs, corresponding to Q — 10.55 GeV and Q — 91 
GeV. In each case, we have allowed for A 5 = 125, 175, 225 MeV. For each of these, we have allowed for 
fi = 2m, m, m/2. 

Notice that the current measurement looks fairly borderline in that is lower than the favored 
range. However, this procedure, which is essentially that of [11]], probably gives too narrow a range of 



predictions, since it is probably not justified to always choose exactly the same value of fio/ m - This 
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Figure 6: Prediction for (^(91GeV)) as a function of the measured (z c (l0.55GeV)), with flb/n^b = 
/i c /m c . A 5 = 125 MeV (dashed line), 175 MeV (solid), and 225 MeV (dashed-dotted). 

procedure underestimates the fi variation, because the values of a one extracts are highly correlated 
with respect to the choice of po/m. So for example, we get a similar prediction if we take /x = rn c to 
extract a and use /io = iTib to determine (zb) as to if we had used /xo = 2m c and 2rrib respectively. 

The optimal choice of jio would correspond to that for which higher order terms are negligible, so 
the subleading order answer, evaluated at the renormalization point //*, is the correct answer. One can 
then ask the question whether r = /i*/m depends on the mass, m. It is straightforward to see that it 
is roughly independent of m if one is performing perturbation theory in a. However, for a leading log 
calculation, r will be formally mass independent only so long as the ratio Q/m\ is much greater than 
mi/ 1712, where Q is the center of mass energy, m\ will be the b quark mass, and m 2 will be that of the 
c quark (assuming common Q). In reality, one probably requires that the scaling between Q and mi 
is much greater than that between mi and 777.2- 

We therefore conservatively allow for independent fio variation in the b and c quark calculations. 
In the case where there are two values of Q involved, we have allowed the full /j,q variation of each. 
In the case where Q = M z for both the b and c quarks, we have always taken both /x greater than 
777 or both less than 777, given that the value of //* is somewhat correlated. This is probably a fairly 
conservative range. The range might be less. 

In Figures 8 and 9, we give the allowed range of possible predictions for (zb) when A5 = 175 MeV, 
given a measurement of (z c ) at Q = 10.55GeV and Q = 91GeV respectively. The corresponding curves 
for A 5 = 125 MeV are slightly more predictive and for A 5 = 225 MeV are somewhat less predictive. 
The allowed range of predictions is greater than in Figure 6 and 7, but still somewhat predictive. 
For example, if we vary (z c ) within 2a around the measured value (^ c (10.55GeV)) = 0.640 we obtain 
0.719 < (* 6 (91GeV)) < 0.782, and similarly for (z c (91GeV)) = 0.487 we get 0.718 < (z 6 (91GeV)) < 
0.763 within 2a. 

Figures 8 and 9 are our main results. They tell us which measurements would be consistent with 
the heavy quark expansion, in light of the perturbative QCD uncertainties. We see that the range is 
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Figure 7: Prediction for (^(91GeV)) as a function of the measured (z c (91GeV)), with //&/m& = /U c /m c 
and the same values of A5 as in Figure 6. 



sufficiently narrow that tests of the heavy quark expansion are possible. A measurement outside this 
range would be an indication of a (A/m) 2 contribution to (z) nonpert (or less likely, an indication of 
nonstandard model physics in low energy QCD scaling). We see that the preliminary numbers do not 
yet require such higher order terms. Reducing the error may or may not leave us inside the range of 
predictions given in Figures 8-9. 



4 Cutoff Moments 

In general, when extracting the moments experimentally, the data does not go down to arbitrarily low 
z. It is therefore necessary to extrapolate the data into this region. This however introduces model 
dependence. The model to which one fits will have some effect on the number which is extracted. An 
uncertainty due to this modeling error can be assigned by attempting to fit to a couple of different 
functions; however this does not necessarily represent the true error. 

In this section, we argue that it might be better to work with truncated moments. By this we 
mean one can define a moment as the integral from zq to 1, rather than from to 1. Ideally, one can 
choose zo as the experimental cut, so long as z is not too large. 

It is however not standard to do this, because unless one integrates to zero, the equation for the 
evolution of the moments does not factorize. In particular, when one defines 



/z 



f 1 dzz N - l f{z), (48) 



and uses 

(49) 

d log u 2 Jz y 
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Figure 8: Prediction for (z&(91GeV)) as a function of the measured (z c (10.55GeV)), for A5 = 175 MeV. 
For each value of the matching scale /i c = 777 c /2 (dashed line), /i c = 777 c (solid line) and /i c = 2?77 c 
(dashed-dotted line) we plot the results for /if, = 7716/2, TTif,, 2/77,6, the higher line corresponding to the larger 

one finds 

dD N r 1 * r 1 

— = / dyy^fiy) / dz^~ x P{z) (50) 

GflOg/^ J Z0 /so /J/ 

In the case z = we see this equation factorizes into the product of moments. However, because of 
the y dependence of the endpoint on the z integration, this is not the case for general zq. The problem 
here is to evaluate the right hand side, you would actually need to know f(y), which is precisely what 
we wish to avoid. What we would like is an expression solely in terms of moments, so we can follow 
the same procedure described in section ||. 

In this section we will focus again on the second moment and we will show that if we make the 
approximation that the evolution equation does factorize, by taking the lower limit on the z integration 
to be Zo, the error we make is actually only of order Zq which is generally rather small. 

So we consider the difference 

pi rl rzo/y 

/ dzzP(z) - / dzzP(z) = / dzzP(z) (51) 

between the function which truly appears in the evolution equation and the approximation one obtains 
by taking the lower endpoint of the integral to be Zq. We consider only the leading log anomalous 
dimension, as the error from the sub leading piece is suppressed by oqcd and should be a small 
correction compared to the error from the leading term. 

As discussed in section |2|, it is sufficient to consider only nonsinglet evolution which involves the 
function 
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Figure 9: Prediction for (,25(91GeV)) as a function of the measured (;Z c (91GeV)), for A5 = 175 MeV. The 
values of fl c and /if, are as in Figure 8, although given that they are somewhat correlated we omit in this figure 
the combinations /i c = T7l c /2, /if, = Iffl^ and jJL c = 2?7l c , /if, = 77lft/2. 



It should be kept in mind that the "+" keeps the function properly normalized. However this should 
be true independent of the arbitrary choice of zq. We therefore take the coefficient of the 5 function 
which needs to be subtracted to be determined from integrating over all z from to 1. (The fact that 
this assumption leads to a smaller error is good confirmation of the expectation that this will better 
approximate the true function). So the error in our approximation is determined from the following 
difference 



where the . . . represents higher order terms in z and we have expanded the logarithm which is obtained 
upon integration. The error is only of order z\. Really, it is of order (zo/y) 2 . However, so long as 
the distribution is weighted towards unity, as it is for a the fragmentation function of a heavy quark 
evaluated at a sufficiently low momentum scale, our statement is approximately correct. 

To check that we do not make too large an error when applying this approximation to a heavy 
quark fragmentation function, we explicitly evaluated the error for different assumed functional forms 
which were claimed to fit data reasonably well in the past. Notice that we are using the function only 
to estimate the error, but not to fit the data. Furthermore, the error estimate only involves those 
values of z which are measured-one never needs the values of the function below zq. So ultimately, it 
is possible to estimate the error based on a function which fits the data, without extrapolation. 

The test functions we used were a Peterson function with a large range of e and also a function of 
the form from ref. namely of the form z a (l — zy . The resulting differences (in percent) between 
the exact and approximate evolution of the second moment are given in Tables 5 (c quark) and 6 (b 




(53) 



which is 




(54) 
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Table 5: Relative difference (in %) between the exact and approximate evolution of the second moment 
for the c quark. 



zo 


0.05 


0.1 


0.15 


0.2 


a = 0.595, (5 = 18.67 




0.15 


0.3 


0.6 


e b = 0.016 


0.05 


0.18 


0.7 


1.3 


e b = 0.049 


0.1 


0.4 


1.0 


2.0 



Table 6: Relative difference (in %) between the exact and approximate evolution of the second moment 
for the b quark. 

quark) for a range of cutoff values, zq. We see that in no case does the error exceed 10% for the c 
quark and 2% for the b quark, and in most of the cases it is considerably smaller. These particular 
numbers correspond to Q — M z , and we have taken A 5 = 200MeV and fi = m, but varying them in 
the ranges we have allowed for in section 2 does not change our conclusions. It is important that the 
error is calculated entirely based on the region above zq, so no assumption is made about extrapolating 
to the unmeasured region. It can evaluated in a fairly model independent manner. 

It is clear that the error is always within the regime of accuracy of our predictions (ie smaller than 
the uncertainty from higher order QCD effects). Therefore, it is as useful to determine the cutoff 
moments as it is to determine the true ones, and probably has a smaller systematic error. We suggest 
the data be presented in the future in this way when there is a cut on z (or as a function of cut 
values for comparison with other experiments). As we have seen, the error when using our procedure 
to extract meaningful numbers should be quite small, of order a few percent. 

5 Conclusions 

In this paper, we have used the measured values of the heavy quark fragmentation second moments 
to test the heavy quark effective theory. By relating c and b quark measurements we have extracted a 
mass suppressed heavy quark parameter. We have also related measurements of c quark fragmentation 
using purely perturbative QCD. The graphs of Figures 1, 8, and 9 are our main results. 

From the purely perturbative analysis, we have seen that the measurements of c quark fragmenta- 
tion at different center of mass energies are consistent with QCD predictions, and tend to favor large 
values of Aqcd- This part of the analysis does not involve heavy quark symmetry and could prove an 
alternative way of determining the strong coupling. 

We have found the value of the nonperturbative parameter a as measured with the c quark lies 
between 79MeV and 422MeV (assuming it is in the overlap of the two measurements) and in the range 
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199MeV to 750MeV for the b quark. A better determination of Aq C d would narrow the allowed range 
of a sizably. For instance, we have obtained that if A 5 = 175 MeV a lies between 208MeV and 396MeV 
for the c quark and in the range 287MeV to 686MeV for the b quark. Although the range is substantial, 
it is of interest to have any measurement at the level of A/m. Furthermore, it is encouraging that 
these values are reasonably small and overlap. With better measurement, one might hope to search 
for deviations at the level of (A/m) 2 . However, given the large perturbative QCD uncertainty, it is 
not yet clear whether this will be possible. 

Finally, we have suggested to use "cutoff" moments, which do not involve extrapolation of the 
experimental data to low z, and therefore will probably have smaller systematic error. We have shown 
that although factorization is not exact, it is a good approximation within the regime of accuracy of 
our calculation. 

There are several measurements which it would be nice to see done or improved. One would want 
CLEO to measure the ratio of D to D* production to test for the importance of higher twist effects. 
One would also want to take advantage of the many c quark states at CLEO to get a better low energy 
measurement of the fragmentation function. A measurement of the production of cccc at LEP would 
be a useful confirmation of the claim that multi c quark production will not distort our predictions. A 
good exclusive measurement on D*'s at LEP would be very useful when using fragmentation functions 
to extract Aq C d- 

It would certainly be advantageous to improve the measurements of heavy quark fragmentation 
functions. Heavy quark fragmentation could prove to be an interesting test of the heavy quark mass 
expansion and a supplement to existing measurements of Qqcd- 
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